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The severity of the sign problem in lattice QCD at nonzero baryon density is measured by the 
average phase of the fermion determinant. Motivated by the equivalence of chiral random matrix 
theory and QCD to leading order in the epsilon regime, we compute the phase of the fermion 
Q_i determinant for general topology in random matrix theory as a function of the quark chemical 

potential and the quark mass. We find that the sign problem becomes milder with increasing 
topological charge. The analytic predictions are verified by detailed numerical random matrix 
fN| simulations. 
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1. Introduction 

The fermion sign problem is a major obstacle for unquenched lattice QCD simulations at 
nonzero density, see [1] for a review at this conference. The severity of the sign problem is 
measured by the average phase of the fermion determinant. Let us denote the Dirac operator by 
D(m;/i), where m and /I are the quark mass and the quark chemical potential, respectively, and 
write detD(m;/x) = re' 9 £ C. Then the average phase of the two-flavor theory is given by (e 2 ' e ). 
It was shown in [2] that the average phase is nonzero for pL < m n /2, whereas for pL > m K /2 it is 
exponentially suppressed in \i 2 F 2 W . Here, m n is the pion mass, V is the volume, and F and £ (see 
(1.1) below) are the familiar low-energy constants in the effective chiral Lagrangian. 

The microscopic regime of QCD (corresponding to the leading order in the £-expansion) is 
defined by the requirement that the scaling variables 



are kept fixed in the V — > °° limit. In this regime exact analytical results for QCD can be derived 
using chiral random matrix theory (RMT). In [3] the average phase was computed for zero topology 
from RMT. An interesting question (e.g., for fixed-topology simulations [4], which could also be 
done at /I ^ using the overlap operator [5]) is whether the sign problem becomes harder or milder 
as a function of the topological charge v. In [6] we have generalized the results of [3] to nonzero 
V. In this contribution, we summarize our main results and present some additional material. 

Before doing so, let us briefly mention important work related to the present subject. The 
leading-order /^-regime corrections to the RMT result for the average phase were computed in [3]. 
The average phase at nonzero temperature (but still at v = 0) was computed from RMT in [7] . The 
distribution of the phase was studied in [8] (see also [9] at this conference), and lattice studies of 
the phase were performed in [10, 1 1], 

2. Random matrix model and sketch of the calculation 

The starting point of our calculation is the non-hermitian random matrix model [12] 



for the Dirac operator, where <I>i and <J>2 are independent complex random matrices of dimension 
(N + v) x N, distributed according to the same Gaussian distribution, P(X) exp(-NtiXX^). The 
matrix in (2.1) has | v| exact zero modes, which allows us to identify v with the topological charge. 
We keep v fixed and take N — > °°. Without loss of generality we assume v > 0; for v < we can 
simply replace v by | v| in our results. 

Going over to an eigenvalue representation of the random matrices <I>i and <3>2, the partition 
function of the random matrix model with Nf flavors becomes [12] 



m = mVL and p. 2 = \i 2 F 2 V 



(1.1) 




(2.1) 



Zy f (a)= / ]Jd 2 z k w v (z k ,zt-,a)\A N ({z 2 })\ 2 l\(m 2 f-zl) 




(2.2) 
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where a = jJ. 2 , A^ is the Vandermonde determinant, and we have introduced the weight function 



4a 



2a 



(2.3) 



that includes a modified Bessel function K. The ensemble average of an observable is given in the 
usual way by including it in the integrand and dividing the resulting integral by Z. 
The phase factor of the squared determinant can be written as 



N W 2 -7 2 
_ TT m h 

det(Dt( jU ) +m ) l = l m 2 



7to_ det(D(ju)+m) 



(2.4) 



with ensemble average 



/ m\ = I det(D(/x)+w) 
\ /w / \det(Dt( jU ) +m ) 



Zy' (a,m) 
Zy f (a) 



Zy f (a)Jc^ 



n^w v (z„4;a)|A w ({ z 2 })| 2 -^^n(^.- z f) 



m- 



' z k /=i 



(2.5) 



iVf+lll* 

where Z v is the RMT partition function with Nf + 1 fermionic quarks and one conjugate 
bosonic quark. The integrals in (2.2) and (2.5) can be done for any N using the orthogonal- 
polynomial formalism developed in [13, 14]. The resulting calculation is too extensive to be re- 
produced here, and we therefore refer to [6] for the details. The main ingredients are orthogonal 
polynomials [12] with respect to the weight function (2.3) and their Cauchy transforms. After tak- 
ing the microscopic limit with the RMT scaling variables m = 2Nm and a = 2Na (corresponding 
to the physical scaling variables in (1.1)) kept fixed, one arrives at 
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(2.6) 



where the subscript s labels the microscopic limit, we have assumed degenerate quark masses, 
Iv,k(z) = z k I v +k{z) with modified Bessel function /, and 
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corresponds to the microscopic limit of the Cauchy transform mentioned above. The calculation of 
this integral is a nontrivial exercise in complex analysis and results in 



1 z ' lx 6A 



4d 



du(-) k u k+l exp 
+ / duu k+l exp 



u \ ( mu\ , , 



8a 



K v 



mu 
4d 



l v+ k{u) 



+A v>k (a,m) (2.8) 



with 



A v h{6c,m) = e 



-2a 



/tov+/t-l i+j< V-l 



V-l-j)!(v+jk-l-7')! 



,f„ !v l / ./)!/!/• 



8d I 



(2.9) 



The v = limit of (2.8) agrees with the result obtained in [3], and A v ,* is a new term that was not 
present for v = 0. This essentially completes the calculation. 

From the above equations one can derive a number of interesting limits, such as the chiral or 
the thermodynamic limit. In the latter limit, we take m and d to infinity keeping their ratio fixed 
and obtain 



2ie th _ j(\-2a/m) Nf+i for2d<m, 
[O for2d>m, 



(2.10) 



where the condition 2d < m corresponds to pL < m K /2. Note that (2.10) is independent of V and 
agrees with the result of [2, 3]. 



3. Results and discussion 



In Fig. 1 we display the dependence of the average phase in the microscopic limit on the 
parameters of the problem. As expected, the sign problem becomes harder as we increase d (i.e., 
jli) and decrease m (i.e., m). 

The dependence of (e^ 1 ® ) on v can be anticipated based on a qualitative argument. At jl = 
the Dirac eigenvalues are purely imaginary, and at \i ^ they move away from the imaginary 
axis into the complex plane. The behavior of the small nonzero Dirac eigenvalues is known from 
RMT [12, 15]. On average, the imaginary part of the small eigenvalues increases with |v|, whereas 
the real part is mainly controlled by jU and depends only weakly on v. Therefore, we expect the 
average phase to increase with | v| . This expectation is confirmed by the analytical results displayed 
in Fig. 1, i.e., the sign problem becomes milder as |v| is increased (as expected, this effect weakens 
if m is increased). This may be an interesting observation for groups doing lattice simulations at 
fixed topology in the e-regime of QCD [16]. As long as the parameter jX in (1.1) is not too large, 
unquenched simulations using reweighting techniques (see Sec. 4 below) are feasible. It would 
then be advantageous to simulate not at v = but at larger values of V. 

As for the dependence on Nf, one would naively expect that more flavors make the sign prob- 
lem harder. However, there is a competing effect, again based on the average behavior of the small 
eigenvalues known from RMT: For small m, Nf has the same effect on them as |v|, i.e., it increases 
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Figure 1: The top row shows the a -dependence of the average phase in the chiral limit for v = 0, 1,2 and 
Nf ■ = 0, 1 , 2. The bottom row shows the dependence on m for fixed a = 1 and the same values of v and Nf. 



the imaginary part but does not significantly affect the real part. As above, this competing effect is 
washed out if rh is increased. The Nf -dependence of (e^' e ) in Fig. 1 is consistent with this qual- 
itative argument. We observe that for v = and small rh, increasing Nf actually makes the sign 
problem milder. Interestingly, for |v| = 1 and in the chiral limit, (e^ e ) is completely independent 
of Nf. For v > the competing effect of Nf does not dominate, i.e., increasing Nf makes the sign 
problem harder as naively expected. 



4. Numerical random matrix simulations 

We have confirmed our analytical results (including finite-./V effects) by extensive numerical 
simulations of random matrices [6]. Examples are shown in Fig. 2 (for the quenched case) and in 
Fig. 3 (for Nf = 2). In the unquenched case, the numerical data were obtained by three different 
reweighting methods: 

1 . reweighting from the quenched ensemble, 

2. reweighting from the phase-quenched ensemble (with r = | detD(m; fx)\ in the measure), 

3. reweighting from the "sign-quenched" ensemble (with r\ cos 6 1 and r\ sin 6 1 in the measure). 

Method 3 minimizes the amount of reweighting that is being done and is therefore expected to lead 
to the smallest statistical errors. This expectation is confirmed in the right plot in Fig. 3. 

It would now be interesting to compare our analytical results with data from lattice simulations. 
We have done initial lattice studies, but our volume was too small to find agreement (more precisely, 
we did not have enough eigenvalues in the microscopic regime). In the future we hope to be able 
to go to larger volumes to confirm the validity of the RMT results for QCD. 
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Figure 2: This figure is for the quenched case (Nf = 0). The top row shows the a-dependence of the average 
phase in the chiral limit for v = 0,1,2. The solid lines are the RMT results for N — >• «>, and the data points 
are from RMT simulations at N = 20. The bottom row shows the dependence on m for fixed & = 1, with 
RMT simulations at Af = 20 and N = 80. The deviations between the analytical results and the data points 
are simply finite-Af effects. The filled area represents the contribution of the A-term specific to v 7^ 0. 




Figure 3: Dependence of the average phase on m at fixed & = 1 (left) and on a at fixed m = 1 (right), both 
for Nf — 2 and v = 1. The solid lines are the analytical results at finite Af [17], and the data points are from 
reweighted RMT simulations at the same value of N. In the left plot sign-quenched reweighting was used. 

5. Summary and outlook 

We have computed analytical results for the average phase of the fermion determinant from 
RMT and found that the sign problem becomes milder with increasing |v|. Our analytical results 
are valid for QCD in the microscopic regime. We have confirmed these results by extensive numer- 
ical simulations of RMT. Interesting projects for the future are the calculation of the distribution of 
the phase for arbitrary topology, lattice simulations at larger volume, and a study of the applicability 
of the sign-quenched algorithm (item 3 above) in other problems. 



6 



Random matrix analysis of the QCD sign problem 



Tilo Wettig 



References 

[1] P. de Forcrand, Simulating QCD at finite density, PoS(LAT2009)010. 

[2] K. Splittorff and J. J. M. Verbaarschot, Phase of the Fermion Determinant at Nonzero Chemical 
Potential, Phys. Rev. Lett. 98 (2007) 031601, [hep-lat/0609076]. 

[3] K. Splittorff and J. J. M. Verbaarschot, The QCD sign problem for small chemical potential, Phys. 
Rev. D75 (2007) 116003, [hep-lat/0702011]. 

[4] JLQCD Collaboration, S. Aoki et. al., Two-flavor QCD simulation with exact chiral symmetry, Phys. 
Rev. D78 (2008) 014508, [arXiv : 0803 . 31 97]. 

[5] J. C. R. Bloch and T. Wettig, Overlap Dirac operator at nonzero chemical potential and random 
matrix theory, Phys. Rev. Lett. 97 (2006) 012003, [hep-lat/0604020]. 

[6] J. C. R. Bloch and T. Wettig, Random matrix analysis of the QCD sign problem for general topology, 
JHEPm (2009) 100, [arXiv: 0812 . 0324]. 

[7] J. Han and M. A. Stephanov, A Random Matrix Study of the QCD Sign Problem, Phys. Rev. D78 
(2008)054507, [arXiv: 0805 . 1939]. 

[8] M. P. Lombardo, K. Splittorff, and J. J. M. Verbaarschot, Distributions of the Phase Angle of the 
Fermion Determinant in QCD, arXiv: 0904 . 2122. 

[9] M. P. Lombardo, K. Splittorff, and J. J. M. Verbaarschot, Fluctuations, correlations and the sign 
problemin QCD, PoS LAT2009 (2009) 171, [arXiv: 0912 . 3109]. 

[10] S. Ejiri, Canonical partition function and finite density phase transition in lattice QCD, Phys. Rev. 
D78 (2008) 074507, [arXiv:0804.3227]. 

[11] J. Danzer, C. Gattringer, L. Liptak, and M. Marinkovic, A study of the sign problem for lattice QCD 
with chemical potential, arXiv: 0907 . 3084. 

[12] J. C. Osborn, Universal results from an alternate random matrix model for QCD with a baryon 
chemical potential, Phys. Rev. Lett. 93 (2004) 222001, [hep-th/04 03131]. 

[13] G. Akemann and A. Pottier, Ratios of characteristic polynomials in complex matrix models, J. Phys. 
A37 (2004) L453-L460, [math-ph / 0404068]. 

[14] M. C. Bergere, Biorthogonal Polynomials for Potentials of two Variables and External Sources at the 
Denominator, hep-th/0404126. 

[15] G. Akemann, J. C. R. Bloch, L. Shifrin, and T. Wettig, Individual complex Dirac eigenvalue 

distributions from random matrix theory and lattice QCD at nonzero chemical potential, Phys. Rev. 
Lett. 100(2008) 032002, [arXiv: 0710 . 2865]. 

[16] JLQCD Collaboration, H. Fukaya et. al., Two-flavor lattice QCD simulation in the epsilon-regime 
with exact chiral symmetry, Phys. Rev. Lett. 98 (2007) 172001, [hep-lat/0702003]. 

[17] J. C. R. Bloch and T. Wettig, to be published. 



7 



